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The stationary solutions of the Gross-Pitaevskii equation can be divided in two classes: those 

> 
-^j- ■ which reduce, in the limit of vanishing nonlinearity, to the eigenfunctions of the associated 

oo ; 

Schrodinger equation and those which do not have linear counterpart. Analytical and numeri- 
cs 

cal results support an existence condition for the solutions of the first class in terms of the ratio 



between their proper frequency and the corresponding linear eigenvalue. For one-dimensional con- 
fined systems, we show that solutions without linear counterpart do exist in presence of a multi-well 
external potential. These solutions, which in the limit of strong nonlinearity have the form of chains 

of dark or bright solitons located near the extrema of the potential, represent macroscopically ex 

O ■ 

cited states of a Bose-Einstein condensate and are in principle experimentally observable. 

•j_J ■ PACS numbers: 03.65. Ge, 03.75.Fi, 47.20.Ky 
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I. INTRODUCTION 

Bose-Einstein condensates of gases of alkali atoms confined in magnetic or optic traps are 
effectively described by the Gross-Pitaevskii equation (GPE), a Schrodinger equation with 
a local cubic nonlinear term which takes into account the interaction among the bosons in a 
mean field approximation [|TJ . The stationary solutions of the GPE represent macroscopically 
excited states of the condensate and have attracted great theoretical interest |2|, [J. The 
existence of some of these states have been also demonstrated in recent experiments. Vortices 
have been observed in two- || or one-component M condensates. Phase engineering optical 
techniques have allowed to generate dark solitons in atomic gases with positive scattering 
length 1,0. 

The excited states observed in 0, ||, |], 0] have linear counterpart, i.e. they are stationary 
solutions of the GPE which reduce, in the limit of vanishing nonlinearity, to the eigenfunc- 
tions of the associated Schrodinger equation ||. However, the GPE can admit also a class of 
stationary states without linear counterpart. These solutions appear for a sufficiently large 
value of the nonlinearity whenever the system is confined in an external multi-well potential 

E- 

In this paper we review the general properties of the stationary solutions of the GPE. In 
Section II we discuss an existence condition for the solutions with linear counterpart in terms 
of the ratio between their proper frequency and the corresponding linear eigenvalue. We also 
show that in the limit of strong nonlinearity these solutions assume the shape of chains of 
dark or bright solitons depending on the repulsive or attractive nature of the interaction. In 



Section |TJ we generalize the asymptotic shape of the states with linear counterpart to find 
a new kind of solutions of the GPE. These correspond to solitons located near the extrema 
of the external potential in a way which generally breaks the symmetry of the system. We 
consider the particular case of a symmetric double-well in Section |TV| and describe all the 
zero-, one-, and two-soliton solutions of this system. In Section we show an example of how 
the solutions without linear counterpart are generated by increasing the nonlinearity. The 
stability properties of the stationary solutions in view of a possible experimental observation 
are pointed out in Section [V]]. 



II. SOLUTIONS WITH LINEAR COUNTERPART 

The stationary solutions of the GPE are defined as 

* li {x,t) = e-l>*M*)> (2- 1 ) 

where \i is called chemical potential, and determined by the equation 

h 2 

- ^ v2 M x ) + U \M X )\ 2 M X ) + V(x)^(x) = n^{x) (2.2) 

with the normalization condition ||Vvll 2 — N\tp^\. For later use we note that the stationary 
solutions are also critical points of the grand-potential functional 
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dx 



EW\-nN[4>\. (2.3) 



Let us consider stationary solutions ip^n which reduce for Uq — > to the eigenfunctions 
4> n of the Schrodinger equation 

h 2 

- —V 2 M X ) + V{x)<j) n {x) = E n (j) n {x). (2.4) 

2m 

This is the linear limit of the GPE which can be obtained, for Uq fixed, by varying the 
norm of the solutions, i.e. the chemical potential jjl. In fact, by substituting in the GPE, 
4>nn{x) = ^jNjjj)x»n( x ), with N n (n) = \\ip^ n \\ 2 and \\x^n\\ 2 = 1, we get 

h 2 

- y~ V 2 x M „(a;) + C/ iVn(^)IX/m(aO| 2 X/m(aO + V(x)x^(x) = yx^x). (2.5) 

For N n (fi) small, the nonlinear term can be neglected and Xim — 4>n- Moreover, we have 

Li~E n + U NM\\(f) 2 n \\ 2 , (2.6) 

i.e. for N n (u.) — > the chemical potential tends to £ n form above or below depending on 
the sign of Uq. 

The above considerations suggest the following existence conjecture. For U > (U < 0), 



solutions with linear limit ip^n ~ J N n (n)(f) n exist only if jj, > S n (jj, < £ n ). Moreover 
N n (fi) — » for jjl —* S n . This conjecture can be i) proved by a theorem in the case n = 0, ii) 
explicitly verified in the solvable case of a one-dimensional system confined in a box, and iii) 
supported by numerical results for multidimensional systems with different potentials M. 



Here, we illustrate point ii) whose results are useful also for the subsequent discussion on 
general external potentials. Let us consider the case of a one-dimensional system confined 
in a box extending from —L/2 to L/2. For Uq > 0, the Jacobi elliptic functions 



tl>^(x) = Am(2(n + l)K(p)(j + ^ 



P , 



where 



K(p) 



-.de 



(2.7) 
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/0 ^Jl-psm 2 9 
is the complete elliptic integral of the first kind with modulus p G [0, 1], and n = 0, 1, 2, . . ., 



solve the GPE under the conditions 
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(2.9) 
(2.10) 



Since K(p) increases monotonously from K(0) = 7r/2, for a given n Eq. ( |2.10| ) has solution 
only if 
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This complies with the conjecture formulated above. 

For Uq < 0, the solutions of the GPE in the box are of the form 

% n (x) = Acn (2(n + l)K(p) (~ + i) + K(p) 

with the conditions 



(2.11) 
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(2.13) 
(2.14) 



mL 2 2 

Since (1 — 2p)K(p) decreases monotonously for p e [0, 1], the n-node solution exists only if 
[i < £ n 'm agreement with the conjecture. 

For both Uq > and C/o < 0, the stationary solutions of the GPE in the box reduce in 
the linear limit to the well known Schrodinger eigenfunctions 
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In the opposite limit of strong nonlinearity, we obtain chains of solitons. For Uq > and 
ji ~^> £ n , we get dark soliton solutions 

VV(X) "^ -p WTT" II tanh \f~ ( x ~ x k)j > ( 2 - 16 ) 

with solitons centered at a^ = (— | + ^ti^) £■ For Uq < and — /x ^> £ n , we have bright 
soliton solutions 
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±(-l) k sech(^^(x-x k )y (2.17) 



with solitons located in x k 
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III. SOLUTIONS WITHOUT LINEAR COUNTERPART 

The soliton chains obtained in Sec. II as stationary solutions of the GPE in a box can 
be generalized in the case of a potential V of arbitrary shape. Let us consider first the case 
Uq > 0. For fi — > oo, the repulsive interaction tends to delocalize the solutions so that a 
Thomas-Fermi approximation holds, i.e. we can neglect the gradient term in the GPE 

U \ij(x)\ 2 ij(x) + V(x)i/>(x) ~ nip{x). (3.1) 

Therefore, the GPE has always the solution 

{ y/{fi-V(x))/U V > V(X) 

[0 /I < V(x) 

and, in the one-dimensional case, also n-node solutions of the form 

ipixn(x) = ip»o(x) II tami ( n ( x ~ x k)j , (3-2) 

provided that the solitons do not overlap, i.e. \xk+i — Xk\ 3> h/^mfi. 

In the attractive case Uq < 0, for /i — > — oo the mean field density tends to be localized 
and the linear potential in the GPE can be neglected with respect to the cubic one 

- ^ 2 i[>(x) + U \iJ>(x)\ 2 i/>(x) ~ ^(x). (3.3) 

In the one-dimensional case, this equation has solutions of the form 



il>v,n(x) = \ JT S Sfc sech * ( x ~ x fc) ' ( 3 - 4 ) 

V ^o fc=0 V ft / 



provided that \xk+i -ii| > h/^—2m/j,. Note that for the corresponding solutions with 
linear counterpart we must have Sk = (— l) fc - However, all sign combinations Sk = ±1 are 
possible in general. 

In order the dark- and bright-soliton chains ( |3.2| ) and ( |3.4j ) to be asymptotic (for fi — > oo 
or fj, — > — oo, respectively) solutions of the GPE, the soliton centers must be extremal points 
of the corresponding grand-potential Q({xk})- Depending on the shape of the external 
potential, stationary solutions without linear counterpart may arise. Consider, for instance, 
the one-soliton solutions. In the strongly nonlinear limit, the soliton width, Tij ^Jm/i for 



Uq > (h/^/— 2miJ, for Uq < 0), vanishes and Q(xx) ~ const — V{x\) (fi(xo) ~ const + V(xq) 
for Uq < 0). If the external potential V(x) is single-well, Q(xi) or Q(xq) will have only one 
extremum. The corresponding solution necessarily reduces in the linear limit to the 1-node, 
for Uq > 0, or 0-node, for Uq < 0, Schrodinger eigenfunction. However, if the external 
potential is multi-well, several one-soliton solutions are possible. Some of them will not 
have linear counterpart, i.e. they disappear when the linear limit is approached. 

The approximate dark- or bright-soliton solutions now introduced make possible a nu- 
merical search for the GPE stationary solutions for any value of the nonlinearity. In fact, 
stationary solutions of a PDE can be found numerically by using a relaxation algorithm 
which converges to the solution "closest" to a given input function. Thus, only solutions of 
which a good approximation is known can be found. In practice, the numerical search of 
stationary solutions can be organized in the following way: 

- choose a sufficiently large |/i| and consider an asymptotic solution of the form of a 
dark- or bright-soliton chain; 

- determine the soliton centers x^ by extremizing the corresponding grand-potential 
tt({x k }); 

- use this approximate solution as input function in the relaxation algorithm; 

- use the obtained output solution as input function in a new relaxation with a smaller 
value of |// 1; 

- repeat the last step until the linear region is reached. 



IV. SYMMETRIC DOUBLE- WELL SYSTEM 

As an example of the general approach outlined in the previous Section, now we determine 
all the zero-, one- and two-solitons solutions for a system confined in a symmetric double- well 
potential. We choose 
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V[x) = m 7 x — muo x + -— (4.1) 



which has two minima in x = ±x m , where x m = Juj 2 /2m'y 4: } and a maximum at x = 0. 

Zero-soliton solutions exist only for the repulsive GPE. For \x sufficiently large, just one 
node-less state is possible. This state extends over the entire double well and is given 
by ipfjjo(x) = J (/J, — V(x)) /Uq for x such that V(x) < [i and ip,jo(x) = otherwise. If 
jjl is smaller than the barrier height uj /A^ A , the above solution still exists and eventually 
approaches the ground state of the Schrodinger equation for the double-well. However, two 
other possibilities appear. We can have a state localized in the left well that joins with 
the identically vanishing solution in the right well and viceversa. These states break the 
symmetry of V(x) and do not have linear counterpart. 



One-soliton solutions are described by Eq. ( |3.2|) with n — 1 in the repulsive case and 
Eq. ( j3.4|) with n = in the attractive one. The corresponding grand-potential becomes a 
function of the soliton centers X\ or xq, respectively. As shown in Fig. 1, for \fi\ sufficiently 
large we have Q(xi) ~ const — V(xi) and fi(xo) ~ const + V(xq). In both cases we have 
three one-soliton solutions corresponding to the three extrema of the external potential. 
The soliton may be found in the maximum or in one of the two minima of the double-well. 
For Uq > 0, the solution of the first kind reduces in the linear limit to the anti-symmetric 
Schrodinger eigenstate with a single node. For Uq < 0, the solution with a bright soliton in 
the maximum of the double-well, even respecting the symmetry of the potential, disappears 
in the linear limit. The two solutions with the soliton, dark or bright, in one of the minima 
±x m of the double- well break the symmetry of V(x) and do not have linear counterpart. 

In the repulsive case, two-soliton solutions are described by Eq. (|3.2| ) with n = 2 and the 
grand-potential becomes the two- variable function Q(xi, £2) whose contour plot is shown in 
Fig. 2. When the distance between the soliton centers is much larger than their width, we 
have Q(xi, 22) — fi(^i) + ^(22). m ^ ne region x\ < X2, fi has a maximum in (— x m , x m ) and 
two saddle points in (0, x m ) and (— x m ,0). We assume that x m ^> Ti/^/mjJ. The stationary 
solution corresponding to the maximum of Q has linear counterpart, namely the symmetric 
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FIG. 1: Grand potential Q as a function of the dark soliton center x\ (left panel) and the bright 
soliton center xq (right panel) for the one-soliton solutions. 
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FIG. 2: Contour plot of f2(xi,X2) for the two-soliton solution in the repulsive case. 



Schrodinger eigenstate with two nodes in the double-well minima. Those corresponding to 
the two saddle points break the symmetry of V(x) and must disappear in the linear limit. 

Other extrema of Q can be found when the centers of the two dark solitons are into the 
same well. As shown in Fig. 2, Q has two maxima in (x m — 5, x m + 5) and (— x m — 5, — x m + 5) 
with 25 > fr/y/mfj,. The corresponding solutions break the symmetry of V(x) and do not 
have linear counterpart. 

In the attractive case, the bright solitons solutions are given by Eq. fl3.4Q . Since the 
GPE is invariant under a global phase change, if we restrict to real solutions for n = 1 (two 
solitons) we have the following two possibilities 






sech (x — xq) J ± sech (x — x\) 



(4.2) 



The contour plots of the functions Q, (xq,x\) obtained by inserting these expressions in 
£l[ij)] are shown in Fig. 3. In analogy with the repulsive case, for both Q + and Q~ we have 
a minimum in (— x m ,x m ) and two saddle points in (0, x m ) and (— x m ,0). The stationary 
states corresponding to the minimum of ft (xo,Xx) have as linear counterpart the lowest- 
energy symmetric and anti-symmetric Schrodinger eigenstates of the double well. Those 
corresponding to the two saddle points break the symmetry of V(x) and do not have linear 
counterpart. 

The functions Q + and Q~ have different behavior when both the soliton centers are inside 
the same well, i.e. for xq ~ x\ ~ ±1. In fact, Q + does not present new extrema while Q~ , in 
analogy with the repulsive case, has two minima in (x m — 5, x m + 5) and (— x m — 5, —x m + 5) 
with 25 •> h/^—2mfi. The corresponding solutions break the symmetry of V(x) and do not 
have linear counterpart. 

An example of stationary states without linear counterpart is shown in Fig. 4. These 
states have been calculated numerically with the procedure outlined in Section |ITJ. The 
parameters used are those of a realistic condensate: m = 3.818 x 10 -26 Kg, uo = 12.75 Hz, 
7 = 10 9 Kg~4m~2s~2 ) U = 1.1087 x 10~ 41 Jm. With these values, the distance between 
the double-well minima is 2x m ~ 92 /xm. 
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FIG. 3: Contour plot of Q + (xq,xi) (left) and Q (xq,xi) (right) for the two-soliton solution in 
the attractive case. The grand potential Q^ is evaluated for the two possible real states 
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FIG. 4: Density of the stationary states ip^,i(x) repulsive and ipti(x) attractive. For comparison, 
we show also the double- well potential V(x). 

V. BIRTH OF STATIONARY SOLUTIONS 



To understand how the GPE solutions without linear counterpart arise by departing from 
the linear limit, consider the following example. We assume a low tunneling regime between 
the two wells of the potential of Section [IV], i.e. ui 3 /h^ A 3> 1, and set 

ip(x) = VN [a (j) (x + x m ) + b (fi (x - x m )] , al + bl = l, (5.1) 

where <fi n {x) is the n-th eigenfunction of the Schrodinger problem with harmonic potential 
\m{2u)) 2 x 2 . Since ip is already normalized to N, for it to be a stationary solution of the 



GPE we have to extremize the energy functional E[ 
small terms, we have 

3 h^ ■ i 



E(b 



Nhu 
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16 u 3 



+ e JwM 1 + 



= £l[ip] + fJ^N [ijj]. Up to exponentially 

3 ^7 4 3u 3 \ 



U" 



2hj 4 



1-68 



NU / m 



2^ V K 



CO 



(l + 2b 4 - 2b 2 



This can be rewritten, up to a constant, as 



E(b ) ~ 6 0V /l-6§ + sign(?7o)— (l + 26, 



tf-26&), 



where 



a; 3 " 3 



h A u 



\ mC/ ( 



(5.2) 



(5.3) 



(5.4) 
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FIG. 5: Energy of the state (|5.l| ) as a function of the parameter bo for different values of N/Nq. 
The upper curve corresponds to the case Uo < 0, the lower one to Uq > 0. For N = the two 
curves coincides. 

The behavior of E(bo) for different values of the ratio N/N is shown in Fig. 5. For 
N <C N , E(b ) has a minimum for 6 = 2^2 and a maximum for b = — 2~2. These extrema 
correspond to the lowest energy symmetric and anti-symmetric linear states. If Uq > 0, for 
iV ~ Nq the maximum at &o — —2^2 bifurcates in a minimum and a maximum. Further 
increasing N, the latter moves toward bo = 0. This describes the birth of the state without 
linear counterpart with dark solitons and its localization into the left well. If Uq < 0, 
the behavior of E(b ) is similar with maxima and minima exchanged. In this case, the 
bifurcation of the minimum at b = 2~? and its move to b = describes the birth of the 
state with 1 bright soliton which localizes into the left well. 

VI. CONCLUSIONS 

We have shown that in presence of an external potential a 1-D GPE can admit stationary 
solutions without linear counterpart. Their existence is strictly connected to the multi-well 
nature of the potential. In the double well example illustrated here, these solutions disappear 
in the limit uo — ► when the potential assumes the shape of a single quartic well. For a 
piece-wise constant double-well, the stationary states here discussed analytically only in the 
limit of strong nonlinearity can be obtained in terms of Jacobi elliptic functions for any 
number of particles in the condensate. 

In H we have also investigated the stability of the stationary states of the GPE under 
different points of view. The results indicate that the soliton-like states, with and without 
linear counterpart, are sufficiently stable on the typical time scales of a BEC experiment. 
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By introducing proper perturbations of the stationary states, a soliton dynamics could also 
be observed. 
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